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Well-posedness of a nonlinear integro-differential problem and its 

rearranged formnlation 

Gonzalo Galiano ^ Emanuele Schiavi § Julian Velasco^ 


Abstract 

We study the existence and uniqueness of solutions of a nonlinear integro-differential 
problem which we reformulate introducing the notion of the decreasing rearrangement of the 
solution. A dimensional reduction of the problem is obtained and a detailed analysis of the 
properties of the solutions of the model is provided. Finally, a fast numerical method is de¬ 
vised and implemented to show the performance of the model when typical image processing 
tasks such as filtering and segmentation are performed. 

Keywords: Integro-differential equation, existence, uniqueness, neighborhood filters, decreas¬ 
ing rearrangement, denoising, segmentation. 


1 Introduction 

This article is devoted to the study of the nonlinear integro-differential problem 

dtu{t,x) = [ )Ch{u{t,y) - u{t,x)){u{t,y) - u{t,x))dy (1) 

Jn 

+ A(uo(x) - u(t,x)), 

u(0,x)=uo(x) (2) 

for (t, x) G Qt = (0, T) X Q. Here, O C (d > 1) denotes an open and bounded set, T > 0, 
A > 0 and uq G BV{Q) nL“(0). The range kernel JCh is given as a rescaling K,h{0 = K,{^/h) of 
a kernel /C satisfying the usual properties of nonnegativity and smoothness. We shall give the 
precise assumptions in Section]^ We shall refer to problem Q-Q as to problem P(0,uo). The 
main results contained in this article are: 

• TheoremThe well-posedness of problem P(n, uq)) the stability property of its solutions 
with respect to the initial datum, and the time invariance of the level set structure of its 
solutions. 

• Theorem]^ The equivalence between solutions of problem P(n, uq) and the one-dimensional 
problem P(r2* *,no*), where n* = (0, |U|), and uq* is the decreasing rearrangement of uq, 
see Section [2] for definitions. 
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• Theorem]^ The asymptotic behavior of the solution of problem P(n*,no*) with respect 
to the window size parameter, h, as a shock filter. 

Problem P(n, uq) is related to some problems arising in Image Analysis, Population Dynam¬ 
ics and other disciplines. The general formulation in Q includes, for example, a time-continuous 
version of the Neighborhood filter (NF) operator: 

NF'"u(x) = -— / e ^ u{y)dy, 

C{x) Jn 

where h is a positive constant, and C'(x) = f^exp (—|u(x) — dy is a normalization 

factor. In terms of the notation introduced for problem P(Q,uo) the NF is recovered setting 
/C(s) = exp(—and A = 0. This well known denoising filter is usually employed in the image 
community through an iterative scheme, 

- u(^Hy))u(^Hy)dy, (3) 

Cn(x) Jn 

withC'„(x) = /Cft (x)-u( (y))dy. It is the simplest particular case of other related filters 
involving nonlocal terms, notably the Yaroslavsky filter [SUES], the Bilateral filter 123 EH], and 
the Nonlocal Means filter [8|. 

These methods have been introduced in the last decades as efficient alternatives to local 
methods such as those expressed in terms of nonlinear diffusion partial differential equations 
(PDF’s), among which the pioneering nonlinear anti-diffusive model of Perona and Malik |19j . 
the theoretical approach of Alvarez et al. [T] and the celebrated ROF model of Rudin et al. |22j . 
We refer the reader to [9| for a review comparing these local and non-local methods. 

Another image processing task encapsulated by problem P(D, uq) is the histogram prescrip¬ 
tion, used for image contrast enhancement: Given an initial image uq, find a companion image 
u such that u and uq share the same level sets structure, and the histogram distribution of u 
is given by a prescribed function T. A widely used choice is T(s) = s, implying that u has a 
uniform histogram distribution. In this case, JC{s) = sign~{s)/s and A is related to the image 
size and its dynamic range, see Sapiro and Caselles |23j for the formulation and analysis of the 
problem. Nonlinear integro-differential of the form 

dtu{t,x)= / {u{t,y) - u{t,x))w{x-y)dy (4) 

Jn 

and other nonlinear variations of it have also been recently used (Andreu et al. [6]) to model 
diffusion processes in Population Dynamics and other areas. More precisely, if u{t, x) is thought 
of as a density at the point x at time t and w{x—y) is thought of as the probability distribution of 
jumping from location y to location x, then u{t, y)w{x — y)dy is the rate at which individuals 
are arriving at position x from all other places and —u{t,x) = — f^u(t,x)w(x — y)dy is the 
rate at which they are leaving location x. In the absence of external or internal sources this 
consideration leads immediately to the fact that the density u satisfies the equation 

These kind of equations are called nonlocal diffusion equations since in them the diffusion of 
the density u at a point x and time t depends not only on u{t, x) but also on the values of u in 
a set determined (and weighted) by the space kernel w. A thoroughfull study of this problem 
may be found in the monograph by Andreu et al. [6|. Observe that in problem P(D,uo), the 


2 






Well-posedness of a nonlinear integro-differential problem and its rearranged formulation 


space kernel is taken as u; = 1, meaning that the influence of nonlocal diffusion is spread to the 
whole domain. 

As noticed by Sapiro and Caselles |23j for the histogram prescription problem, and later 
by Kindermann et al. [I6] for the iterative Neighborhood filter or by Andreu et al. [6] 
for continuous time problems like Q, these formulations may be deduced from variational 
considerations. For instance, in m, the authors consider, for u £ the functional 

Jiu)= g{u(x) - u{y))w{yL-y)dxdy, (5) 

JflxU 

with an appropriate spatial kernel w, and a differentiable hlter function g. Then, the authors 
formally deduce the equation for the critical points of J. These critical points coincide with the 
fixed points of the nonlocal filters they study. For instance, if g{s) = )Ch{\/a)da and rc = 1, 
the critical points satisfy 

= 7^ I ^hiu{x) - u{y))u{y)dy, 

Jfi 

which can be solved through a fixed point iteration mimicking the iterative Neighborhood filter 
scheme Q. On the other hand, choosing g{s) = s (or some suitable nonlinear variant) and 
considering a gradient descent method to approximate the stationary solution, equation (Q is 
deduced. Similarly, g{s) = |s| and w = 1 leads to the histogram prescription problem. 

Although the functional ([^ is not convex in general, Kindermann et al. prove that when K, 
is the Gaussian kernel then the addition to J of a convex fidelity term, e.g. 

J{u-Uq) = J{u) + A||U - Uo\\\2(^y 

gives, for A > 0 large enough, a convex functional J, see |T6l Theorem 3.1]. 

Thus, the functional J may be seen as the starting point for the deduction of problem 
P(n,tto), representing the continuous gradient descent formulation of the minimization problem 
modeling Gaussian image denoising. Notice that although the convexity of J is only ensured for 
A large enough, the results obtained in this article are independent of such value, and only the 
usual non-negativity condition on A is assumed. 

The outline of the article is as follows. In Section]^ we introduce some basic notation and the 
definition of decreasing rearrangement of a function. This is later used to show the equivalence 
between the general problem P(n,Mo) and the reformulation P(fl*,Mo*) in terms of a problem 
with a identical structure but defined in a one-dimensional space domain. This technique was 
already used in [12] for dealing with the time-discrete version of problem P(fl,uo)j in the form 
of the iterative scheme Q. See also |13lll4j for the problem with non-uniform spatial kernel. In 
Section]^ we state our main results. Then, in Section]^ we introduce a discretization scheme for 
the efficient approximation of solutions of problem P(n,uo), and demonstrate its performance 
with some examples. In Section]^ we provide the proofs of our results, and finally, in Section]^ 
we give our conclusions. 

2 The decreasing rearrangement 

Given an open and bounded (measurable) set C {d > 1) let us denote by jnj its Lebesgue 
measure and set fl* = (0, |II|). For a Lebesgue measurable function tt : fl —)■ M, the function 
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g G M —)■ mu{q) = |{x G : m(x) > (?}| is called the distribution function corresponding to u. 
Function niu is, by definition, non-increasing and therefore admits a unique generalized inverse, 
called its decreasing rearrangement. This inverse takes the usual pointwise meaning when the 
function u has not flat regions, i.e. when |{x G n : u{x.) = (?}| =0 for any q € In general, 
the decreasing rearrangement u* : fl* —)• M is given by: 

{ ess sup{rt(x) : x G 0} if s = 0, 
infjg G M : mu{q) < s} if s G fl*, 

essinf{ri(x) : x G 0} if s = |n|. 


Notice that since u* is non-increasing in fl*, it is continuous but at most a countable subset of 
fi*. In particular, it is right-continuous for all a G (0, |n|]. 

The notion of rearrangement of a function is classical and was introduced by Hardy, Little- 
wood and Polya [15] . Applications include the study of isoperimetric and variational inequalities 
ISnilZlIITKIB], comparison of solutions of partial differential equations [SHlEllMlIinillllll], and 
others. We refer the reader to the monograph [21] for a extensive research on this topic. 

Two of the most remarkable properties of the decreasing rearrangement are the equi-measurability 
property 


/ fiu{y))dy 

Jn 



f{u^{s))ds, 


( 6 ) 


for any Borel function / : M —>■ ]R_|_, and the contractivity 


u* - 'y*||LP(r2*) < ||w - 'y||LP(n)) 


(7) 


for u,v £ LP{Q), p G [l,oo]. 

For the extension of the decreasing rearrangement to families of functions depending on a 
parameter, e.g. t G [0,T], we first consider, for t fixed, the function u{t) : H —)• M given by 
u{t){x) = u{t,x), for any x G H. Then we define w* : (0,T) x —>■ M by s) = n(t)*(s). 

3 Main results 

Our first result ensures the well-posedness of problem P(n,Mo) for initial data with 

bounded total variation. In addition, we show that the level sets structure of the solution is 
time invariant. Before stating our results, we collect here the main assumptions on the data 
problem, to which we shall refer to as (H): 

• 0 C is an open, bounded, and connected set {d> 1). 

• The final time, T, which simulate the time horizon of the diffusion process is a real, positive 
fixed number. 

• The parameter A is a real, nonnegative fixed number. 

• /C G is nonnegative. 

• tto G A := n BV{^1) is assumed to be, without loss of generality, non-negative. 
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Basic facts but also advanced results about the space of bounded variation BV{Q) can be found 
in the book by Ambrosio et al. [5]. Notice that, depending on the space dimension d > 2 we 
have the continuous injections BV{Q) When d = 1 we have X = BV{^). 

Theorem 1 Assume (H). Then there exists a unique solution u G C°°([0, T]; A) of problem 
P{Q,uo). In addition, ifuoi,uo 2 £ X and ui,U 2 G C°°{[0,T]] X) are the corresponding solutions 
to problems P{Q,uoi), P{Il,uo 2 ) then 

ll^^l - 'W2||L°°(0,T;L2(f7)) < C'll'^Ol “ U02\\L^{n), (8) 


for some constant C > 0. 

Finally, suppose that uo(xi) = uo(x 2 ) for some xi,X 2 G II. Then u{t,xi) = u{t,X 2 ) for all 

t G (o,r]. 

Remark 1 The existence and stability results of Theorem may be extended to more general 
zero-order terms in the equation 0 of problem P{Q,uo). For instance, we can consider a 
function f : [0, rjxOxM —satisfying /(•, x, s) G L°°(0, T), f{t,-,s) G BV (fi), and f{t,x,-) G 
This regularity coincides with the initially obtained for the integral term of equation 
0 in the approximation procedure to construct the solution. In addition, i/uo(xi) = uo(x 2 ) 
implies f{-,xi,-) = /(•,X2,-), then the time invariance of level sets holds. 

Replacing the set H by 0* and the initial data uq G T by uq G A* = BV{fl^), Theorem 
ensures the existence of a solution of problem P(n*,uo*)- Observe that fl* C M is bounded 
because 0 C is bounded (assumption (H)) and this implies BV{It^) C and X^ = 

BV{n^). 

In the following result we obtain some properties of solutions of the one-dimensional problem 
P(n*,uo). Although the corollary is valid for any interval in M, we keep the notation fl* for 
simplicity. The corresponding result for the discrete-time version, with A = 0, of problem 
P(n*,no*) may be found in [12]. 

Corollary 1 Assume (H), and let v G C°°{[D,T]-,X^) he the solution of problem P(n*,uo), for 
some nonincreasing vq G T* . Then 

1. sigii{dsv{t, •)) = sign(9sUo) < 0 a.e. in 11*, for all t G (0,T). 

2. For t > 0, v{t, 0) < uo(0) and v{t, |n|) > uo(|l^|)- 

3. If JC is odd then v{t, s)ds = vo(s)ds, for t > 0. 

4. IfJC G andvo G , for m G N andl <p < oo, thenv G C^{[0,T]-,W^’P{n^)). 

5. // A = 0 and JC G C^(M) is such that -|- K,h{f) > 0 then v{t, •) —>• const, as t ^ oo. 

Remark 2 Condition in point 3 is a natural symmetry condition for convolution kernels and 
it is satisfied, for instance, by the Gaussian kernel. Condition in point 5 is also satisfied by the 
Gaussian kernel, if h is large enough. 

The next result establishes the connection between problems P(ll,uo) and P(n*,uo*). 
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Theorem 2 Assume (H). Then, u G T]; V) is a solution of P{TI,uq) if and only if 

u* G C'“([0, T]; V*) is a solution o/P(n*,no*). 

Theorem implies that the solution of the multi-dimensional problem P(n,uo) may be 
constructed by solving the one-dimensional problem P(n*,uo*)- Indeed, using the level sets 
invariance asserted in Theorem we deduce 

u{t, x) = u^{t, s) for a.e. x G {y G 0 ; uo(y) = tto*('S)}; 


for all t G [0,r]. When image processing applications are considered, by property 1 of Gorol- 
lary HI the solution to P(n,Mo) may be understood as a contrast change of the initial image, 

Uq. 

Indeed, this property also implies that if, initially, uq has no flat regions, and therefore uq* is 
decreasing, then the solution of P(n*, uq*) verifies this property for all t > 0. Then, Theorem 
implies that the solution of P(n,uo) has no flat regions for all t > 0. 

The last theorem is an extension of a result given in m for the discrete-time formulation with 
A = 0. In it, we deduce the asymptotic behavior of the solution u* of problem P(Q*,uo*) (and 
thus of u of problem P(Q,uo)) in terms of the window size parameter, h. Although we state it 
for the Gaussian kernel, more general choices are possible, see na Remark 2], 


Theorem 3 Assume (H) with lC{ff) = e~^^ and uq G T having no flat regions. Suppose, in 
addition, that uq* G Then, for all {t,s) G [0,r] x 12*, there exist positive constants 

01,02 independent of h such that the solution tt* G C'°°([0, T]; o/P(n*, uq*) satisfies 


dtufit,s) = A(uo*(s) - s)) -h a]kh{t,s)h^ - 02 


(9) 


with 


and with oi 


kh{t, s) 


lCh{u^{t, |n|) - u^{t,s)) 

\dsU^{t, |n|)| 


l/(2-y/T), and 02 ~ I. 


lCh{u^{t,Q) - u^{t,s)) 
|clsu*(t,0)| 


( 10 ) 


Two interesting effects captured by Q are the following; 

1. The border effect (range shrinking). Function kh is active only when s is close to the 
boundaries, s ^ 0 and s ^ |fl|. For s 0, kh{t,s) < 0 contributes to the decrease of the 
largest values of u* while for s ~ |n| we have kfi{t, s) > 0, increasing the smallest values of 
tt*. Therefore, this term tends to flatten tt*. In image processing terms, a loss of contrast 
is induced. 


2. The term 

dssU*{t,s) 

\dsU^{t, s)p 

is anti-diffusive, inducing la^e gradients on u*(t, •) in a neighborhood of inflexion points. 
In this sense, the scheme (IM is related to the shock filter introduced by Alvarez and 
Mazorra [2] 

Vt T F{GcrVxxi G(jVx^Vx — 0, (II) 
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where is a smoothing kernel and function F satisfies F{p,q)pq > 0 for any p,q € 
Indeed, neglecting the fidelity, the border and the lower order terms, and defining F(p, q) = 
we render ® to the form (11). 


p 

qIqI 


This property can be exploited to produce a partition of the image so the model can be 
interpreted as a tool for fast segmentation and classification. An example is proposed in 
the numerical experiments where a time-continuous version of the NF is implemented. 


4 Discretization and numerical examples 

For the discretization of problem P{Q,uo), for uq : C —)• M, we take advantage of the 

equivalence result stated in Theorem Thus, we first calculate a numerical approximation, 
uo*, to the decreasing rearrangement uq* : C M —)• M and consider the problem P (O*, uq^ ). 

Then, we discretize this one-dimensional problem and compute a numerical approximation, 
v:[0,T] X O* —>■ M. By Theorem]^ v is, in fact, an approximation to u*, where u : [0, T] x O —)• M 
is a solution to problem P(n, uq). Then, we finally recover an approximation, u, to u by defining 

u{t, x) = v{t, s) for a.e. x G {y G Q : 'Uo(y) = 'Uo*('S)}. (12) 

Inspired by the image processing application of problem P(ll,uo), we consider a piecewise 
constant approximation to its solutions. Let 11 C be, for simplicity, a rectangle domain 
and consider a uniform mesh on H enclosing square elements (pixels), Tmn, of unit area, with 
barycenters denoted by x^m for m = 1,..., M and n = 1,..., A^. Given uq G L°°(11) n BV (11), 
we consider its piecewise constant interpolator no(x) = Mo(xmn) if x G Tmn- 

The interpolator uq has a finite number, Q G N, of quantized levels that we denote by qi, 
with max('Uo) = Qi > ■ ■ ■ > QQ = min(Mo). That is uo(x) = Ylf=i QjXEji^), where Ej are the 
level sets of ho, 

Flj = {x G H : ho(x) = qj}, j = 

Since ho is piecewise constant, the decreasing rearrangement of ho is piecewise constant too, and 
given by 

Q 

uo*{s) = '^qjXij{s), (13) 

1=1 

with Ij = [aj-i,aj) for j = and oq = 0, oi = \Ei\, 02 = |Lli| \E 2 \,...,aQ = 

\E,\ = |Si|. 

Let u be a candidate to solve problem P(H*,ho*). Due to the time-invariance of the level 
sets structure of the solution to this problem, see Theorem we may express v as 

Q 

= ^Cj{t)xijis), (14) 

1=1 

with ci(l) > ... > CQ{t), for t G (0,T], Cj(0) = = qj, for j = Substituting v in 

equation ([^, we get, for s G Ij and j = 1,... ,Q, 

Q 

= '^^h{ck{t) - Cj{t)){ck{t) - Cj{t))p,k + Mcj - Cj(t)), (15) 

k=l 
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with /Xfc = Ofc — afc_i = \Ek\. Since, by assumptions (H), the right hand side of (15) is Lipschitz 


continuous, the existence and uniqueness of a smooth c = (ci,. .., cq) : [0, T] —>■ satisfying 

(15) and c(0) = follows. 

For the time discretization, we take a uniform mesh of the interval [0,r] of size r > 0, and 
use the notation = c(t„), with tn = nr, and n = 0,1, 2,... Then, we consider the following 
implicit time discretization of problem (15). For j = 1,..., Q and n > 1, solve 


Q 


cj = cr‘ + r ^ Kn(cl - - C])n + tA(c“ - c“). 


(16) 


k=l 


Since problem (16) is a nonlinear algebraic system of equations, we use a fixed point argument 
to approximate its solution, c"', at each discrete time tn, from the previous approximation 
Let Then, for m > 1 the problem is to find c"’’™ solving the linear system 


Q 


n,m n— 

S- = S' 




n,m—1 


— C 


n,m—Iw n,m 


K n.m n,7n\ . \ / U n,m\ 

Cfc “S + 'rA(c,-- C.’ ), 


n,m\ 

rj ’ 


(17) 


k=l 


for j = 1,... ,Q. We choose the stopping criterion ||c”’™ — c^’^^^Hoo < tol, for values of tol 
chosen empirically, and then set c” = c"'’™. 

Finally, using formula (12), the expression of the initial datum (13), and the definition (14), 
we recover a piecewise constant approximation to the original problem, P(n,rto), taking 

u{t,x) = dj if t £ [tn,tn+i), X £ {y G n : Uoiy) = Qj}. 


4.1 Example. Histogram based image segmentation 

As an application we consider a Grand Challenge in Biomedical Image Analysis. This is a 
computer vision problem in biomedicine which consists of overlapping cells segmentation and 
subcellular nucleus and cytoplasm detection, see |26], [25]. The dataset was downloaded from 
the Overlapping Cervical Cytology Image Segmentation Challeng^ ISBI 2014. 

The data set is composed by 512 x 512 real and synthetic images containing two or more 
cells with different degrees of overlapping, contrast, and texture. The phantom images allow the 
quantitative analysis of segmentation procedures through their ground-truth, which is carried 
out by using the Dice similarity coefficient, DC: for two sets (images) A and B, 


DC = 


2|AnB| 

\a\ + \b\- 


Observe that values of DC close to one indicate high coincidence of the images, that is, of the 
ground-truth segmentation and the segmentation obtained with our method. 

For running our algorithm, that is, providing an approximation, c"", of (15), we consider 


the usual number of image quantization levels, Q = 256. The fidelity term is ignored (A = 0), 
and the range window parameter, h, is set as h = 25 for nucleus detection, and as h = 5 for 
cytoplasm detection. The tolerance in the fixed point loop ([ItI) is taken as tol = l.e — 5. As a 


'http: //cs . adelaide . edu.au/~carneiro/isbil4_challenge/index.html 
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Table 1: Example: Segmentation results for some samples of testOO dataset. 


Sample 

1 

15 

30 

45 

60 

75 

90 

All (mean) 

Cytoplasm DC 

0.99 

0.99 

0.99 

0.99 

0.99 

0.99 

0.99 

0.98 

Nucleus DC 

0.93 

0.94 

0.82 

0.88 

0.90 

0.81 

0.85 

0.87 

Execution time 

2.44 

4.33 

5.11 

5.12 

5.57 

5.94 

4.18 

5.33 


stopping criterion, we consider a combination of a maximum number of time iterations (1000), 
and an energy stabilization criterion, 

|J(c")-J(c”-i)| < l.e-10, 

where J{c^) is the discrete version of the functional given by (§, for u) = 1 and g{s) = 
exp(— s//i). Finally, we implement a variable time step, T(n), inspired by the proof of exis¬ 
tence of solutions and given by, for re > 2, 


I J(c” — J(c"' 2)1 

with r(l) = r(0) = (— c^)gk}) ^■ In the experiments, we observed that r(re) 
ranges from order 10“^ in the first iterations to order 10“^ just before convergence. 

We summarize our results for the testOO dataset in Table where we show the DC for some 
specific samples, and the mean DC of the ninety samples contained in the dataset. We may 
check that DC values are very high for the segmentation of both regions of interest (cytoplasm 
and nucleus), always above the range obtained in [26L . The execution times are given for a 

Matlab implementation of the algorithm, running on a standard laptop (Intel Core 17-2.80 GHz 
processor, 8GB RAM). 

In Figure we show the segmentation process for the two regions of interest. The first 
column corresponds to the initial image. The second column, to the background extraction, and 
the third column to the nucleus segmentation. Thus, the cytoplasm is the difference between the 
images shown in the third and second column. Finally, the fourth column shows the difference 
between the ground-truth nucleus segmentation and the obtained with our method. 


5 Proofs 

Proof of Theorem We divide the proof in several steps. 

Step 1. Existence of a local in time solution to an auxiliary problem with smooth data. 

We assume uq G VF^’“(f2), and consider the following auxiliary problem, obtained using the 
change of unknown u = rce^* in ([^, for some positive constant /r to be fixed: 

dtw{t,x) = [ fCh{e>^^iw{t,y)-w{t,x))){w{t,y)-w{t,x))dy 
Jn 

+ X{wo{x) - w{t,x)) - gw{t,x), (18) 

for (t,x) G (0, To) X n, and for the initial data t(;(0, ■) = wq = uq £ 1T^’°°(11). Here, Tq > 0 will 
be fixed later. 
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Time discretization. Let V G N, r = Tq/N and tj = jr, for j = 0,...,N. Assume that 
Wj G is given and consider the functional A : L°°(n) —)• L°°{Q) given by 

^ (wj{^) + T {wj{y) - Wj{^))){ip{y) - ip{^))dy 

+ tXwo{x)^, 


for X G 0. Observe that if A has a fixed point ip, then we may define rcj+i = ip to get the 
following semi-implicit version of (18) 


«^j+i(x) =Wj{yi) + T )Ch{e^^^^{wj{y) - Wj{^))){wj+i{y) - Wj+i{^))dy 

Jn 

+ tX{wo{x) - Wj+l{x.)) - TpWj+lix.). (19) 


We have, 


|A((y9(x)) - A(V^(x))| = 


< 


1 + r/i 

2r|ll| 


JCh{e>^^^iwj{y) - Wj{x))) 

X |<^(y) - V’(y) - (<^(x) - V'(x))|dy 

\\)Ch\\oo\\T - V'lloo- 


1 + r/i 

Therefore, for p > 2|n|||/C/i||oo, the mapping A is contractive in L°°(0*), and a unique fixed 
point, Wj+i verifying (19) does exist. 

following ur 

ll^i'j+llloo < 


We have the following uniform estimates for tCj+i. One one hand, from (19) we obtain 

1 


l + '7‘(A + /r — 2|0| 11/C/ill 
which gives (recall p > 2|0|||/C/i||oo) the uniform estimate 

||'(«i+l||oo < Mo 


^ 11 Wj 11 oo + 'rA||u/o|| OO^ 5 


( 20 ) 


with Mo depending only on ||t(;o||oo- 

On the other hand, since wo,Wj G 1T^’°°(0), we deduce from (19) rcj+i G 1T^’°°(0). This 
regularity allows to differentiate in ( |19[ ) with respect to the fc—th component of x, denoted by 
Xk, to obtain for a.e. x G H, 


Fi(x 


dw 


'i+i 


dxi 


, , ^ , X dwi 


. X , X 


with 


Fi(x) =l-br(^A + ^-b f ICh{e^^^{wj{y) - Wj{x 


F 2 (x) =1 - re^*^ / IC'h{e>^*^{wj{y) - Wj{x))){wj+i{y) - Wj+i{x))dy, 
Jn 
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from where we deduce 

^ 1 + r{\ + ^,) + (^ + 2re^'^ |f^|Mo||/C;||oo) ||Vu;,-||oo) • 


Solving this differences inequality, we find that, by redefining ^ to satisfy fj, > 2e^'^°\Q\ ||/C'^||ooMo, 
we obtain the uniform estimate HVrCj+iHoo < Mi, with Mi depending only on UVu^oHoo- This 
election of fj. is possible by restricting To to be 


To < 


1 


log 


2|0|||/C;^||ooMo’ 


( 21 ) 


Time interpolators and passing to the limit r —)> 0. We define, for (f,x) G {tj,tj^i] x Q, the 
piecewise constant and piecewise linear interpolators 


W^'^\t,X.) =Wj+l{-X.), = Wj + l{x.) + — \wj{-x.) - Wj + l{x.)). 

T 

Using the uniform L°° estimates of Wj^i and Vwj^i, we deduce the corresponding uniform 
estimates for ||and implying the existence 

of to G -L°°(0, To; VU^’°°(n)) and w G W^’°°{Qtq) such that, at least in a subsequence (not 
relabeled), as r —>■ 0, 

^ w weakly* in L'^(0, To; 1U^’~(U)), 

^ w weakly* in (22) 


In particular, by compactness 

—)> w uniformly in C'([0,To] x U). 


Since, for t G (tj,tj+i], 

|u;(^)(t,x) -'u}(^)(t,x)| = —^(u;j(x) -u;j+i(x))| < T\\dtw^'^^\\L<->{QT^), 

we deduce both w = w and 


—)> w uniformly in C'([0,To] x Q). 


(23) 


Considering the shift operator arW^'^\t, ■) = wj, and introducing the approximation Ot* = , 

for t G {tj,tj+i], we may rewrite (19) as 


{t, ^)= ICh (e(l*(u;(") {t, y) - (t, x))) (u;M {t, y) - (t, x))dy 

Jn 

+ A(uio(x) - w^'^\t,:s.)) - p.w^'^\t,x.), 


(24) 


and due to the convergence properties (22) and (23), we may pass to the limit r —)• 0 in (24) to 
deduce that u; is a solution of (18). 


Continuation of the solution to an arbitrary time T. Given the solution, w, of problem (18) in 
Qto, we may consider the same problem for the initial datum w{Tq, •). Since w{Tq, •) G kU^^(U) 
and the constant Tq > 0 only depends on |n|, ||/C)j||oo and ||mo||oo, see (20) and (21), we obtain a 


11 















G. Galiano, E. Schiavi, and J. Velasco 


new solution w G C{Tq, 2Tq] Clearly, this procedure may be extended to an arbitrarily 

fixed T. Once this is done, a boot-strap argument allows us to deduce w G (^“(O, T; lV^’“(n)), 
implying that u = rce^* G (^“([O, T]; is a solution of P{Q,uo) in Qt- 

Step 2. Non smooth initial data. 

Let us consider a sequence uqs G such that, as e —)■ 0, 


uoe —>• uq in L°°{Q,) (25) 

TV(uo), (26) 


where TV denotes total variation with respect to the x variable. Let us denote by to the 
corresponding solution of P(0,uoe)- 

First, notice that is uniformly bounded in L°°{Qt) with respect to e as a consequence of 
estimate (20) and property ( [2^ . We then obtain directly from equation 0 that 


dtUs is uniformly bounded in L°°{Qt)- 


(27) 


Since uoe is smooth, we may deduce an bound for Vug as in Step 1, not necessarily uniform 
in e, but which allows us to differentiate equation Q with respect to Xk- After integration in 
(0, t), we obtain 


= ^^^{x)Geit,x)(l +X f {Ge{T,x)) ^dx), (28) 

OXk OXk ^ Jo ' 

with Gir(t, x) = exp ( — (A -|- r]^(T, x))dT ), and 

%(Lx) = {K:'i,(ue(t,y) - Ue(t,x))(ue{t,y) - Ueit,x)) (29) 

JCh{us(t,y) - Ueit,x))'^dy. 


Since Kh G 1V^’“(M), we have % uniformly bounded in L°°(Qt) and so and G^ Therefore, 
using (26) we deduce from (28) that 


Vug is uniformly bounded in L°°{0,T; L^(Q)). 


(30) 


Bounds (27) and (30) allow to deduce, using the compactness result [zH Cor. 4, p. 85], the 
existence of u G C([0,r]; V) such that Ug —)• u strongly in LP(Qt), for all p < oo, and a.e. in 
Qt- 


Similarly to the smooth case, this convergence allows to pass to the limit e —)• 0 in ([^ 
(with u replaced by Ug) and identify the limit tt as a solution of P(Xl,uo). Again, the property 
u G L°°(Qt) and a boot-strap argument leads to u G G°°(0,T; A). 

Stability and uniqueness. Let uoi) ^02 £ BV(Pl) and ui, U2 G G°°([0, Tj; A) be the correspond¬ 
ing solutions to problems P(n,ttoi); P(Pl,uo2)- Set u = ui — U2 and uq = uiq — U2o- Then u 
satisfies 


dtu{t,x) = j (^^(ui{t,y) - ui(t,x)) - 4>(u2(t,y) - U 2 (t,x))'^dy 
+ A(uo(x) - u(Lx)), 

u(0,x) =no(x). 
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for (t,x) G Qt, with <l>(s) = ICh{s)s. Multiplying this equation by u, integrating in n and using 
the Lipschitz continuity of (with constant Cl) and Young’s inequality, we deduce 


dt I \u{t,x)\'^dx<CL I I \u{t,y) - u{t,x)\u{t,x)dydx+ ^ f |Mo(x)|^dx 


n Jn 


A 


|u(t,x)pdx 


<Cl(^I |M(t,x)|(ix^ +{Cl-^) I \u{t,x)\^dx 


+ ^^ko(x)pdx. 


Finally, using Jensen’s and Gronwall’s inequalities, we deduce •)IIl 2 (q) < C'||no||L 2 (f^) for 
all t G (0, T), and the result follows. 

Time invariance of the level sets. The proof of this property is similar to the proof of the 
stability property. Let uq G Y and u G T]; Y) be the corresponding solution to problem 

P{Q,uo). Assume uo(xi) = uo(x 2 ), and set Ui{t) = u{t,Xi), i = 1,2. Then, from equation 0 
we get 


dt{ui{t) - U2it)) = j (^^{u{t,y) - ui{t)) - ^>(u(t,y) - U2it))'^dy 
- X{ui{t)) - U2{t)). 


Then, the Lipschitz continuity of <h and Gronwall’s lemma allow us to deduce the result. □ 


Proof of Corollary To prove point 1, notice that from (in dimension d = 1) we deduce 


sign(9sUe) = sign(uQg) < 0 a.e. in (0,T) x fl*. 


a property that also holds in the limit e —)■ 0. Point 2 of the theorem follows from evaluating 
equation © in s = 0 and s = |n|, using that v{t, •) is decreasing for all t > 0, and Gronwall’s 
inequality. Point 3 is a consequence of the assumption on the symmetry of /C, under which the 
integral term in vanishes when it is integrated in U*. Point 4 is easily deduced by successive 
derivation of (1^) (which also holds for e = 0, under regularity assumptions). Point 5 is again 
deduced from (|28]) and the decreasing character of and v. Since, TV{ve{t,-)) —)• TV{v{t, ■)) 
and, using point 2, TV{v{t, •)) < c for all t > 0, we have that the integral term in is evaluated 
inside a closed interval. Therefore, using the assumptions of point 5, we get %(t,s) > C 2 > 0 
uniformly in (t, s). Finally, we obtain the result from (28) in the limit t —)• oo and e —?■ 0. □ 


Proof of Theorem We split the proof in two steps. 

Step 1. First we treat the case in which uq has no flat regions, that is when |{x G : u(x) = 
g}! = 0 for any g G M. By the invariance of the level sets structure proven in Theorem we 
deduce that neither the solution u of P{Q,uo) has flat regions. Then mu{t,-) and u^{t,-) are 
strictly decreasing, implying u^{t, mu{t, q)) = q for any q G M. According to [211 Theorem 9.2.1], 
we have Sju* = dg^ where 

f Sn 

ip{t,s)= / — (Lx)dx, (31) 

■J{uit)>u4t,s)} Ot 
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and we used the notation {u{t) > = {y G : u(t,y) > u^{t,s)}. Integrating Q in 

{u{t) > u^{t, s)} we get 


(/?(t,s) = 


^h{u{t,y) - u{t,x)){u{t,y) - u{t,x))dydx 


' {u{t)>Ut{t,s)} Jo, 

+ A / no(x)dx-A / u{t,x)dx = Ii+I 2 +h- 

J {u{t)>Ut{t,s)} J {u(t)>Ui,{t,s)} 

Due to the u and u* level sets equi-measure, it is immediate that 

I 3 = —X / n*(t, a)da. 

Jo 

The equi-measurability property implies 


(32) 


(33) 


h = 


J{u{t)>Uf{t,s)} Jflt 

from where we deduce 


/ /Cft(u*(t, cr) — M(t, x))(n*(t, cr) — u{t,x))dadx, 

J r2* 


h = 



JCh{u^{t, a) — u^{t, r))(u*(t, a) — u^{t, T))dadT. 


(34) 


0 JO* 


To deal with the term I 2 we observe that due to the invariance of the level set structure, as 
stated in Theorem we have that, for all t G [0,T] and s G D*, there exists a G D* such that 

{x G D : u{t,x) > tt*(t, s)} = {x G D : Uo{x) > uo*(a)}. 

Recalling that u and uq have not flat regions and taking the measure of these sets we deduce 
s = a. Therefore, 


l 2 = X uo{x)dx= / uo^{a)da. 

J {uo>uot,{s)} Jo 


(35) 


Finally, substituting in identity (32) the expressions (31), (33), (34) and (35), and differentiating 
with respect to s, we deduce the result. 

Conversely, let u be a solution of P(D*,uo*). Since uq has not flat regions, Uq^ < 0 in D*, 
and by point 1 of Corollary we have dsv{t,s) < 0 in [0,T] x D*. We define 

u{t,x) = v{t, s) for a.e. X G L(s) = {y G D : uo(y) = uo*(s)}, (36) 

and for all t G [0,T]. Observe that since uq has not flat regions, we have |T(s)| = 0 for all s. 
Therefore, since dgV < 0, we also deduce that u has not flat regions. By construction, 

|{x G D : u{t,x) > u(t,s)}| = |{x G D : Uo(x) > Uo*('S)}| = s, 


implying u* = u. Differentiating in ( |36[ ) with respect to t and using that u is a solution of 
P(D*,mo *)5 we get, for x G L{s), 

dtu{t, x) = dtv{t, s)= }Ch{v{t, a) - v{t, s)){v{t, a) - v{t, s))da + A(uo(s) - v{t, s)) 

Jn, 

= / ICh{u^:{t,a) — u{t,x)){u^:{t,a) — u{t,x))da + X{uo{x) — u{t,x)) 

JUt 


= / }Ch{u{t,y) - u{t,x)){u{t,y) - u{t,x))dy + X{uo{x) - u{t,x)), 

Jn 
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where we have used again the equi-measurability property Q. 

Step 2. We now treat the general case in which uo € X may have flat regions. We use the 
following lemma. 

Lemma 1 Let uq G X. Then there exists a sequence uqj G X such that uoj has no flat regions 
and uoj uq in X. 

We may then apply the Step 1 of this proof to each uqj to obtain that Uj G (^“([O, T]-,X) is 
a solution of P(n,uoi) (without flat regions) if and only if {uj)^ G C°°(0,T;d4) is a solution of 
(uoj)*)- Now we perform the limit j —)> oo. 

Let u G C°°{[0,T]] X) and v G C'°°(0,T;df*) be the solutions of problems P{Ll,uo) and 
P(n*,uo*) ensured by Theorem 

Using the strong continuity of the decreasing rearrangement operation in see 0, and 

the stability property ([^ applied to problem P{Ll,uo), we obtain 

||m* - (Wj)*||Loo(0,T;L2(n*)) < \\u “ IIi,°°(0,T;L2(q)) < UiHuo - Uoj\\L2(n)- 

The same arguments in reverse order applied to problem P(n*,uo*) leads to 

\\v - {Uj)*\\L°°(0,T-,L2{n^)) < C'2||uo* - {uoj)*\\L2(n^) < <^211^0 - U0j\\L2{n)- 

Therefore, using the triangle inequality we deduce 

||u - t(*||L°o(0,r;L2(n^)) < (Cl + C2)||uo - MOj ||l,2(fi) —)■ 0, 

as j —>■ oo. □ 

Proof of Lemma [T1 

In this proof, we rename uq by u and uoj by Uj. Let, foi i G I , Ei = {x G Ll : u(x) = Qi} with 
\Ei\ > 0, be the collection of flat regions of u which is, at most, countable. Thus, I C N. Let 
Xb. and P{Ei) denote the characteristic function of the set Ei and its perimeter, respectively. 
We consider the functions 

j , s ^ Xg,(x) min(gi - Qj+i, 1) 

“ P{jil+P{Ei))+v{x))' 

where v G BV(LI) is a non-negative function without flat regions. Observe that since u G BV (fl) 
we have P{Ei) < oo for all i G E 

Consider, for j G N, the sequence of L°°{Q) functions 

. r u{-k) if X e Q\ Ujg/ Ei 

n, (x) = u(x) - ^ y,](x) = <^ ^ ^ ^ 

i&I I ^‘^{j{^+P{Ei))+v{x)) 

We have: 

(1) Uj has no flat regions in O. Let g G M. We use the decomposition 

{x G Ll : Uj{x) = q} = {x G Ll\ Uj Ei : u{x) = g} Ujjx G Ei : Uj{x) = g}. 

If g = gi for some i G I then x G Ei, and, by definition, Uj{x) = q^ if 

1 

+ P{Ei)) + v{x)) 
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which is not possible. Therefore, ii q = qi we have \uj = qi\ =0. 11 q ^ qi for all i G / then 
|{x G n\ Uj Ei : u{x) = O'} I = 0, so 

\uj{x) =q\ = |r;(x) = -j'(l + T’(Ei)) + min(gi - qi+i, l)/{i^iqi - g))| = 0, 
since v has no flat regions. 

(2) Uj u in L'P{Q) for any 1 < p < oo. This is immediate, since |«(x) — rtj(x)| < i. 

(3) Uj G BV{Q) and Uj u in BV{Q). According to [5l Proposition 3.38], for each i G I, we 

can find a sequence G with 0 < < 1 such that —)> Xe^ in as h —)• 0, and 


hm [ \Vwi\ = TV{xe,) = P{E^) < oo, 

since u G BVi^L). We also introduce a regularizing sequence Vh G (7°° (12) such that Vh > 0 and 
Vh^vin BV{n). Let gl = wlmm{qi - qi+i,l)/{i^{j{l + P{Ei)) + Vh)). Then, 



1 

Pj{l + P{Ei)) 


[ |V<| + 

Jn 


1 

{ij{l + P{EiW 



(37) 


implying that is uniformly bounded in BV{^}) with respect to h. Therefore, there exists 
g'^ G BV{Q) and a a subsequence of (not relabeled) such that 5 )^ —)■ 5 * strongly in L^(n) as 
/i — 7 - 0. Since, by the Dominated Gonvergence Theorem we have g^^ —>• in L^(12) as /i —?• 0, we 
deduce g^ = pj G BV{Q). Taking the limit /i —> 0 in (37) we get 


^ I + P{EiW In 


(38) 


with c > 0 independent of i and j. Thus, using the definition of uj, the triangle inequality, and 
(38) we get 

TViuj) < TV{u) + Y,TV{p>^ < TV{u) + 

iei ^ 

Therefore, Uj G BV{ld). Finally, from the definition of uj and (38) we have 

TV {u — Uj) < TV (pi) —)■ 0 as j —>■ 00 . 
i&I 


□ 


Proof of Theorem [S]. Define 


I{t,s)= / IChiu^{t,a) — u^{t,s)){u^{t,a) — u^{t,s))da. 
JUt 


(39) 


Since uq has not flat regions, we have Uq^ < 0. Then, due to points 1 and 2 of Corollary we 
have dgU^ < 0 in [0,T) x D*, and u*(t,D*) C tto*(12*) for all t G [0,r], respectively. 

Let us consider the inverse of rt*(t, •), the distribution function of u, m„(t, •). Using the 


change of variable s = mu{-, z) and writing a = mu{-, q), we obtain from (39) 

rut{t,0) 


rn*(t,0) 

Ii{t,z) := I{t,mu{t,z)) = / IChiq - z){q - z) 

Ju^(t.\Q\) 


dq 


\dsU^{t,muit,q))\ ■ 


(40) 
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Using the explicit form of JC and integrating by parts, we obtain 


h{t, z) = — [ kh{mu{t, z)) + / IChiq - z) 


{dsU^{t,muit,q)))^ 


dq^, 


(41) 


with kh given by (10) 


By assumption, function 


fit,q) = 


dssU*{t,mu{t,q)) 


{dsU^it,muit,q)))^ 

is bounded in [u^{t, (|U|)), u*(t, 0)] and by point 4 of Corollary it is continuously differentiable 
in (u*(t, |U|),0)). 

Consider the interval = {q : \q — z\ < Vh}. By well known properties of the Gaussian 
kernel, we have 


and 


K{h) := / K,h{q-z)dq< / K,h{q)dq = hy/i:, 

J Jh 

JCh{z — q) < e“^/^ if q ^ Jh = {q ■ \q ~ ^ Vh}. 


In particular, from (43) we get 


'jF 


K^h{q - z)f{t,q)dq 


< 0{h°') for any a > 0. 


(42) 

(43) 

(44) 


Taylor’s formula implies 
fu 4 o) 


^h{q- z)f{t,q)dq= / IChiq-z){f{t,z) + 0{\fh))dq 

n*(|Q|) Jjh 

+ [ K.h{q- z)f{t,q)dq. 


Therefore, from (41), (44) and (42) we deduce, using < 0, 

/r . .. ,, d‘lsU*{t,mu{t, z)) 

\dsU^{t,mu{t, z))\‘^ 


h{t, z) = —{ khimuit, z)) - 


K.{h) + 


Then, the result follows from (40) substituting z by u*(t,s). □ 


6 Conclusions 

In this paper we studied a general class of nonlinear integro-differential operators with important 
imaging applications, such as the denoising-segmentation Neighborhood filtering. 

Although the corresponding PDE problem is multi-dimensional, we showed that it can be 
reformulated as a one-dimensional problem by means of the notion and properties of the decreas¬ 
ing rearrangement function. We proved the well-posedness of the problem and some stability 
properties of the solution, as well as the equivalence between the multi-dimensional and the 
one-dimensional solutions to the problem. 
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Some other interesting properties were deduced for the rearranged one-dimensional version of 
the problem, such as the time invariance of the level sets of the solution (inherited by the multi¬ 
dimensional equivalent solution), and the asymptotic behavior of the solution as a shock-type 
filter. 

Future work will point to the use of rearranging techniques for the generalization of the 
model to include nonlocal effects induced by non-homogeneous spatial kernels, like in equation 
S- As already showed for the discrete time problem m, this situation is much more involved 
suggesting the consideration of the relative rearrangement functional. 
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Figure 1: Example. Nucleus and citoplasm segmentation process. First column corresponds 
to the initial image. Second column, to the background extraction, and third column to the 
nucleus segmentation. The citoplasm is the difference between the images shown in the third and 
second columns. Finally, fourth column shows the difference between the ground-truth nucleus 
segmentation and the obtained with our method. 
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